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Cosmological gamma ray bursts (GRBs) are the brightest explosions in the Uni- 
^ | verse. Satellite detectors, such as Beppo-SAX, HETE2 and more recently Swift, 

have provided a wealth of data, including the localization and redshifts of subsets 
of GRBs. The redshift distribution has been utilized in several studies in attempts 
to constrain the evolving star formation rate and to probe GRB rate evolution 



in the high-redshift Universe. These studies find that the GRB luminosity func- 
tion and/or the rate density evolve with redshift. We present a short review of the 
problems of constraining GRB rate evolution in the context of the complex mix of 
biases inherent in the redshift measurements. To disentangle GRB rate evolution 
from the biases prevalent in the redshift distribution will require accounting for 
the incompleteness of the observed redshift sample. We highlight the importance of 
formulating a 'complete GRB selection function' to account for the main sources of 
bias. 
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1 Introduction 



Cosmological gamma ray bursts (GRBs) are the brightest explosions in the 
Universe. They show a bi-modal distribution in d urations of short bursts o f 



19931 ) 



around 0.3 s length and long bursts of around 30 s (jKouveliotou et al. 
Their ffuence (time-integrated flux) and time- variability implies that their in- 
ner engines must be compact and highly energetic. The association of some 
long-duration or 'classical' GRBs with supernovae implies that they are asso- 
ciated with the core collapse of massive stars. The few short bursts that have 
been well localized show no supernova association and are associated with both 
early and late type galaxies — see Berger (2006). Short bursts most likely have 
a different evolutionary channel compared to long GRBs, one possible scenario 
being the merger of neutron star-neutron star or neutron star-black hole bi- 
naries. In the following discussion, we focus on observations of long GRBs, 
hereafter just GRBs, from the suite of satellite detectors, and the redshifts 
obtained from a subset of localized GRBs. 



Because of their extreme luminosities in 7-rays, GRBs are a promising new 
probe of the high-z Universe. They provide the opportunity to study their 
local interstellar environment, host galaxies and the intergalactic medium. The 
association of GRBs with short-lived massive stars implies that they may track 
the evolving star formation rate (SFR). These exciting possibilities provide 
strong motivation to use the measured GRB redshift distribution as a tool to 
constrain the SFR and GRB rate evolution. 

To constrain GRB rate evolution requires an accurate model for the GRB lumi- 
nosity function (LF), the selection function of the GRB satellite detectors and 
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the influence of using different ground-based telescopes and techniques with 
varying quality to measure redshifts from the GRB afterglow. Complicating 
the analysis further is the assumption that GRBs represent one unified popu- 
lation with a distribution of luminosities encompassed by a single luminosity 
function. Determining the effects of these biases on the redshift distribution 
is fundamental to constraining the astrophysical GRB parameters. Here, we 
provide an overview of recent work on constraining GRB rate evolution and 
the SFR using the measured GRB redshift distribution. We describe from 
first principles the flux-limited GRB selection function and follow this with a 
discussion on the 'complete' GRB selection function, which encompasses all 
selection effects but is presently poorly understood. 



1.1 The evolving star formation rate 



Several observation-based estimates of the evolving SFR have been made in 
recent years, mostly using the observed rest-frame UV luminosity density of 
the galaxy population. However, there has not been a consensus amongst 
workers on the form of SFR evolution because of systematic uncertainties, 
such as the amount of dust extinction, and because of the difficulties of the 
observations. Ear ly studies of the UV and IR l uminosity densities o f star 



forming galaxies (ILilfy et al. 



1996 



Madau. Pozzetti fc Dickenson 



Madau et al 



1996 



Connolly et all 11997 



19981 ) suggested that the SFR inc reases steeply 



Steidel et al 



from z = to z ~ 1 — 2, and then declines at higher z. However. 
(119991 ) favour a near-constant SFR density at z ~ 1 — 4, and model a mono 



tonically increasing rate out to z — 10. Contrary to some of these results, 
Sullivan et al.l (120001 ) find a more gentle increase in the rate. 
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Hopkins fc Beacoml (120061 ) have constrained the SFR by fitting to ultraviolet 
and far- infrared data obtained up to 2006. When these data are taken together, 
they show that it is possible to constrain the SFR out to z ~ 6, with especially 
tight constraints for z < 1. Figure 1, taken from their work, shows that the 
SFR is well constrained from z = — 1 but is not well defined at z > 2, as 
shown by the larger error bars and smaller density of data. These data indicate 
that the SFR may be either roughly constant or declining at z > 5. 



Because of the large uncertainty in the SFR at z > 1, simple 



parametrized 



mode 



s wit h 3 free parameters 



(120001 ) and 



r ave b een used as fitting functions. 



Porciani fc Madau 



Madau fc Pozzetti 



(1200 ll ) provide observation-based parametrized 
SFR density models for three scenarios. The first matches most measured 
UV continuum and Ha luminosity densities and includes an upward correc- 
tion for dust extinction. This model increases rapidly in z — — 1, peaks in 
z = 1 — 2, and gently declines at higher redshifts. A second model remains 
approximately constant at z > 2 and the third model is based on studies 
that suggest that the evolv i ng SF R up to z ps 1 may have been overestimated 



(jCowie. Songaila fc Barger 



19991 ). while the rates at high z m ay have been 



sever ely underestimated due to large amounts of dust extinction (IBlain et. al. 



199J). 



The model functions of 



Madau fc Pozzettil (I2000h and 



Porciani fc Madaul (120011 ) 



have been supplemented by more recent fits to the data (cf. Figure 1 ). Eq ua- 



tion ([T]) shows the functional form employed by 
their fit to the SFR: 



Hopkins fc Beacoml (120061 ) in 



where the evolution function e(z) is defined so that e(0) = 1. 



The observed GRB redshift distribution can potentially be used to constrain 
the independently derived observation-based SFR density, assuming that the 
GRB rate directly follows the SFR. Selection effects, such as dust obscura- 
tion in the host galaxies, wh ich affect the observed SF R, may also modify the 



observed GRB distribution. iBlain fc Natarajanl (120001 ) show that the ratio of 



GRBs with and without optical transients provides a n important new p r obe o f 



Djorgovski et. al 



(120031 ) 



the star formation history of the Universe. Similarly, 
discuss the likely association of observed GRB afterglows with dust obscu- 
ration in the host galaxies. In the following section we review the data and 
techniques employed to probe the link between the GRB and SFR distribu- 
tions. 



2 Analysis of the GRB redshift distribution 



2. 1 Pre-Swift and Swift redshifts 

From 1997 to 2007 February, about 487 GRBs have been localized using satel- 
lite detectors linked to rapid follow-up ground-based optical telescopes, be- 
ginning with Beppo-SAX in 1997 and the High Energy Transient Experiment 
(HETE). The Swift satellite, a multi- wavelength GRB observatory, has local- 
ized about 186 GRBs up to 2006 Nov and a further 35 have been localized by 
HETE and the satellite detector INTEGRAL. 

From the 487 localizations, 99 redshifts have been obtained with reasonable 
confidence and 4 of these are associated with short GRBs. Several 'anomalous' 
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GRBs — namely GRBs 980425, 030329, 031203 and 060218— were observed at 
relatively small redshift (0.0085, 0.168, 0.105, 0.0331 respectively) and were 
also under-luminous compared to the mean luminosity of long GRBs. Among 
these four, only the extremely under-luminous GRB 980425 and very long 
duration GRB 060218 are really peculiar, while the other two may belong 
to the faint end of the burst distribution. It has also been suggested that 
GRBs 030329 and 031203 could be accounted for if they are normal bursts 
seen off-axis, or they may belong to a sub-population of under-energetic GRBs 



( 1 Coward! 



2005|). 



Figure 2 plots the cumulative redshift distribution of 42 and 48 pre- Swift and 
Swift bursts respectively. It is clear that the distributions are biased differently. 
As of 2006 November, the median redshifts of the two samples are about 1.0 
and 2.4 respectively. The discrepancy has been partly explained by comparing 
the trigger criteria of the BAT detector onboard Swift with BATSE. Band 
(2006) found that the longer temporal accumulation for the trigger of the 
Swift BAT makes it more sensitive to long-duration bursts. Also the BAT 
detects bursts in a lower-energy band than BATSE did, and long-duration 
bursts are softer, on average, than short-duration bursts. Hence, the BAT 
preferentially detects long-duration bursts; because the burst durations are 
time dilated from cosmic redshift, the BAT is more sensitive to high-z GRBs 
than was BATSE. Further to this, accurate localizations provided by Swift 
and faster response times of telescopes with spectroscopic capabilities have 
led to more efficient redshift determinations. 
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2.2 The GRB luminosity function 



Band (2006) points out that the detectors, BAT and BATSE, select two differ- 
ent subsets of the GRB population. Unifying the two distributions is fraught 
with difficulty and relies on a complete knowledge of the detector selection 
function and an accurate model for the GRB luminosity function. Unfortu- 
nately the GRB luminosity function is uncertain. Recent studies using the 



Swift redshift distribution have attempted to constrain both 
tion and GRB rate evolution. For example, iFirmani et al 



uminosity func- 



()200J) applied a 



joint fit to 33 GRB spectroscopic redshifts and 220 redshifts inferred from 



an e mpirical variability-luminosity relationship (IFenimore fc Ramirez-Ruiz 



20001 ) with the aim of constraining both GRB luminosity and the SFR as func- 
tions of z. They found that the best fit to these data implies that the luminosity 
is best modelled using a double power-law with slopes —1.6 and —2.6. They 
also find that the best fit to the data is achieved usin g an approximately linear 



evolut ion of the GRB luminosity. In another analysis, 



Guetta. Piran fc Waxman 



( 120051 ) derived an empirical double power-law luminosity function assuming 



parametrized SFR models, but obtained slopes of —0.1 and —2. 



An ex ample LF is shown below using the double-power law form from lGuetta. Piran fc Waxman 



(120051 ). They use lower and upper luminosity limits, 1/Ai and A 2 , respectively, 
to define the isotropic (not corrected for beaming) luminosity function, /(Lj SO ), 
of GRB peak luminosities L iso , in the interval logL iso to logL iso + d\ogL iso : 



/(Aso) = < 



(L iso /L*) a , L*/A x < L iso < L* 



(L iso /L*)?, L* < L iso < A 2 L* , 



(2) 
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where / f(Li SO )dLi SO is normalized to unity. Following lGuetta. Piran fc Waxman 



( 120051 ). we employ their parameter values a = —0.1 and (3 = —2 for the LF 
slopes, with A x = 30 and A 2 = 10 and an isotropic-equivalent break peak 
luminosity of L* « 3.2 x 10 51 erg s _1 . We use this luminosity function in the 
next section to calculate the GRB scaling function. 

There are many alternative LFs that have been employed to constrain the 
slope of the LF and the SFR at high z. For example Daigne et al. (2006) 
use a single power law for their Monte Carlo simulations of the GRB redshift 
distribution, and obtained different slopes assuming 3 different SFR models 
and GRB peak energy distributions. Assuming the peak GRB energies are 
described by an Amati relation^ or by a log-normal distribution, they find a 
power law slope of —1.5 to —1.7 with minimum and maximum luminosities of 
(0.8 — 3) x 10 50 and (3 — 5) x 10 53 erg s -1 respectively. They point out that 
the observed luminosity distribution is biased towards brighter events because 
the detectors are not sensitive enough to detect the low-luminosity GRBs. 



3 The GRB flux-limited selection function 

The selection function is a fundamental feature of any flux-limited survey with 
instrumental bias. It incorporates knowledge of the distribution of luminosities 
of the 'target sources' and the sensitivity of the detector. An accurate model 
of the flux-limited selection function allows one to separate out intrinsic as- 
trophysical parameters, given enough statistics. This is well known in optical 
1 A correlation between the isotropic equivalent GRB energy and the corresponding 
peak frequency in the source rest frame — see Amati et al. (2002); Atteia et al. (2003); 
Lamb et al. (2004). 
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and radio surveys of large numbers of galaxies and is equally important in 
relation to the GRB data obtained from Swift and other satellites if the data 
are to be used for cosmological studies. 



The fraction of all bursts potentially detectable by a satellite depends on many 
factors, including detector area, detector efficiency, angular field of view and 
instrumental noise. In the absence of instrumental selection effects, such as 
the temporal sensitivity to GRBs (e.g., Band 2006), the detection fraction 



depe nds completely on the f 



e.g., ISchaefer. Deng fc Band! (EOOlh ; 



ux sensitivity thresho l d of th e instr ument-see 



Schmidt 



(|2Q0l|); 



Norrisl (120021 ) 



3. 1 The Band function 

The GRB photon spectral luminosity can be approximated by the Band et al. 
(1993) parametrized fit, expressed as: 



S{E) = A x 



E 



100 keV 



exp 



(P-a)E 



E h 



E <E b 



( E b 



(3) 



( E 



exp ((3 -a), E> E b 



V VlOOkeV/ VlOOkeV, 
with S(E) the source rest-frame photon spectral (or differential) luminosity, so 
that A is in photons s -1 keV -1 ; S(E)dE is the number of photons emitted per 
second in the energy interval E to E + dE. Integrating ES(E)dE over some 
energy band then gives the total (all angles) photon bolometric luminosity 
in that band. Assuming isotropic emission, dividing by the surface area of a 
sphere centered on the source converts this to the photon flux in that band, 
with units photons s -1 cm"" 2 . The energy spec tral indices, a and j3 can be 
obtained from the BATSE spectral catalogue of . 



Preece et al. 



200(1 



(120001 ) and have 
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typic al values —1 and —2.25 respectively, with a spectral break at Ef, = 511 
keV ( INatarajan et. al.l . 120051 ) but we note that any value for Ef, ~ 500 keV is 
a reasonable approximation. 



3.2 Flux limited GRB detection 



As shown by lGuetta. Piran fc Waxmanl (120051 ). the flux limit, -Pum^um, z), for 
a GRB of peak luminosity L lim (z) and photon spectral luminosity S(E) at a 
bolometric luminosity distance d^z), is 



L lim S$$%S(E)dE 
^d 2 L (z) J* S(E)dE 



-Plim(-^-'lim) Z) 



(4) 



where E\ < E < E 2 defines the detected energy band and the (1 + z) factor in 
the integration limits accounts for the redshift of the rest frame source spec- 
trum. Typical values for the flux sensitivity limit, Pu m , of the detectors BATSE 
and Swift are 10~ 7 and 10~ 8 ergs cm -2 s _1 respectively with corresponding 
sensitivity bands 50 — 300 and 15 — 150 keV. 

The luminosity distance used above is related to the proper distance dp(z) by 
d L (z) = (l + z)d P (z) , (5) 

and the proper distance is defined by 

2 

H dp(z)/c— Jdx/h(x, Qyi) ■ (6) 
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where H is Hubble's constant and the dimensionless Hubble parameter h(z, Qm) 
is given by 

h(z, n M ) = ^i-n M + (i + z ) 3 n M (7) 



for a spatially flat cosmology. 

The fraction of detectable GRBs, or scaling function (those observed with 
peak flux > Pu m ), is 



where L lim (Pi im , z) is obtained by solving equations (jlj). The GRB scaling 
function, Eq. ([8]), is the inverse of the detector selection function. Figure 3 
plots iPgkb{. z ) using the LF defined by equation ([2]) with slopes a = —0.1 and 
(3 = —2 and using the Swift detector flux limit 10 -8 ergs cm -2 s -1 . The curve 
highlights how GRB detection evolves strongly with redshift. 

The GRB scaling function is fundamental to any model used to fit to GRB 
redshift data. It encompasses both the GRB luminosity function uncertainty 
and the detector response to a spatial distribution of transient fluxes. But it 
does not include other selection biases from either the satellite detectors or 
the redshift measurements from ground-based telescopes. This issue will be 
discussed further in Section 14.21 

3.3 The GRB probability distribution function 

There are a variety of approaches in modeling and simulating the GRB redshift 
distribution. A straightforward method that highlights the key parameters is 



oo 




'ISO ) 
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described below. The all-sky event rate distribution of transient cosmological 
events, such as core-collapse supernovae and GRBs, can be expressed by the 
event-rate equation (e.g. Coward, Burman & Blair 2001): 



dR/dz = ATT(c 3 r /H 3 )e(z)F(z)/(l + z) 



(9) 



where R(z) is the all-sky event rate, as observed in our local frame, for sources 
out to redshift z. The factor (c 3 ro/-£/o 3 ) has the dimensions of inverse time, as 
does the event rate R; the factor r is the "local" rate density of events, the 
event rate per unit volume at the present epoch. The (1 + z) denominator in 
([H]) accounts for the time dilation of the observed rate by cosmic expansion, 
converting a source-count equation to a rate equation. Equation flU]) represents 
the event-rate in shells of redshift dz; its form is defined by the expression 



e(z)F(z)/(l + z). 

The function F(z) is determined by the cosmological model. For flat-A cos- 
mologies (Peebles 1993, p. 332): 



We assume a flat-A cosmology (Qm = 0.3, = 0.7) with H = 70 km s 



The probability distribution function for the potentially observable GRB rate 
variation with z can be expressed as 



F(z) = [H d P (z)/c] 2 /h(z,nu) ■ 



(10) 



Mpc 



-i 



p(z) = ^ GRB (z)e(z)F(z)/(l + z) , 



(11) 
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with normalizing factor 

^max 

N= J iJ GRB (z)[e(z)F(z)/(l + z)]dz . (12) 
o 

Figure H] plots equation (TTTT) using the same scaling functions as shown in 
figure For definiteness we assume GRB rate evolution traces a SFR model 
that peaks in z = 1 — 2, and remains approximately constant at z > 2. 
Equation pip with the normalizing factor N' 1 replaced by 47r(c 3 ro/-£^o 3 ) e 
is the differential GRB rate as observable by a satellite with a flux limited 
scaling function V'grb; the e factor defines the field of view and efficiency of 
the detector and e(z) is the evolving GRB rate density normalized to unity at 
z = 0. 



3.4 Utilizing the GRB redshift distribution 



The fact that GRB redshifts have been measured from the local Universe 
to z = 6.9 provides strong motivation for using the redshift distribution to 
constrain the evolution of the GRB rate density Previous studies have focussed 
on co nstraining GRB rate evolution by assuming that they trace the evolving 



SFR. 



Firmani et al. 



(12004 ) find that their best fit for the SFR, from a joint 
fit using 33 GRB spectroscopic redshifts and 220 redshifts inferred from an 
empirical variability- luminosity law, implies a general shape that is consistent 
with other independent tracers of the SFR, such as the rest-frame luminosity 
densities of galaxies. 

The study by Daigne et al. (2006) employed the log N - logP and peak energy 
distribution of BATSE GRBs and the HETE2 fraction of X-ray rich flashes 
to constrain GRB evolution. Instead of fitting directly to these distributions, 
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Daigne et al. used Monte Carlo simulations to study the effect of the uncer- 
tainties discussed above on the GRB rate. They find a GRB to supernova 
rate ratio of 1:10 5 — 10 6 , assuming no beaming factor, but do not attempt 
to quantify the uncertainty in this range. The most significant result of their 
analysis is the prediction of a positive GRB evolution that exceeds the best 
fit models for the evolving SFR. They quote a GRB rate density at z = 7 of 
6-7 times larger than at z = 2. Daigne et al. interpret this as evidence for 
GRBs not tracing the evolving SFR but rather tracing the evolution of GRB 
progenitor rate efficiency This is an attractive idea given that several studies 
indicate that some GRBs are associated with dwarf irregular low-metallicity 



galaxies (IFruchter et al. 



20061 ) 



Another study, by iTruong fc Dermerl (120061 ). uses a GRB model with a non- 



evolving luminosity function to fit to the redshift and jet opening angle distri- 
butions to constrain GRB evolution. The authors use seven variable param- 
eters in their data fits, with four of them partly constrained by observation. 
Their best fit to a luminosity function, assuming a single power law, finds a 
slope —3.25, an intrinsic beaming factor of 34 — 42 and a GRB rate of one 
in every 600,000 years in a Milky Way type galaxy. Similarly to Daigne et 
al., they find positive GRB evolution at high z, but their calculated GRB to 
supernovae ratio of about 10 4 is 1 — 2 orders of magnitude smaller. Truong & 
Dermer point out that they could only find statistically acceptable joint fits 
to the data if the GRB rate density is increasing monotonically at z > 5, but 
they do not attempt to provide confidence limits to the fits. 



14 



4 Open issues 



4-1 An evolving GRB rate density 1 ', 



Daigne et al. (2006) and iTruong &: Dermerl (120061 ) reach similar conclusions 
for the evolution of the GRB rate density; they find that the o bserved GRB 



Firmani et al. 



(120041 ) reach a sim- 



rate at high z is higher than that expected, 
ilar conclusion, but explain this higher than expected rate as evidence for an 
evolving GRB luminosity function. These results are very interesting and may 
be important to our understanding of how GRBs formed during earlier cosmo- 
logical epochs. If they are robust, there are several possible reasons why the 
GRB rate is relatively excessive at high- z. The number density of GRB pro- 
genitors may have been enhanced in earlier cosmological epochs because of the 
expected higher density of massive stars in low-metallicity host environments. 



Very rec ent work has shown th at long GRBs are more probable in low-metallicity 
galaxies: 



Fruchter et al 



(120061 ) used Hubble Space Telescope images of 42 GRB 
afterglows and their host galaxies in z < 1.2, and found that dwarf irregular 
galaxies are the predominant host galaxy type. In addition they found that 
the GRBs in their sample are located in the brightest parts of their hosts, not 



where the b 



Stanek et al 



ue lig 



it characteristic of active star formation is most intense. 
(120061 ) have shown that host galaxies of the five long GRBs with 
redshifts z < 0.25 have distinctly low oxygen abundance in comparison with 
star-forming galaxies in 0.005 < z < 0.2 from the Sloan Digital Sky Survey. 



Dwarf irregular galaxies undergo prolonged but gradual star fo rmation, and 
are predicted to contain less than 10% of their metallicity in stars (jCalura &: Matteucci 
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20061 ) . For dwarf irregulars, the stellar metal density is a slowly increasing 
function throughout their histories, implying a continuous but relatively low 
SFR. If GRBs are more probable in low-metal enriched environments, GRB 
evolution should track the evolving low-metallicity component of the SFR. 

There is also the strong possibility that the GRB population consists of several 
intrinsically different populations. The observation of the so-called 'anomalous' 
under-luminous GRBs, namely GRBs 980425, 030329, 031203 and 060218, 
may imply a more numerous population based on their relatively small red- 
shifts. Does this imply a different class of progenitor from that of the classical 
GRBs? What are the implications for the low-energy tail of the luminosity 
function? 



4-2 The complete GRB selection function 

Knowledge of the GRB selection function, ?Pgrb( z ), described by Equation 
(jHJ), is defined by a GRB LF, a flux-limited sample and a choice of cosmol- 
ogy dependent luminosity distance. It would represent the 'complete selection 
function' in the absence of all other instrumental selection biases, which un- 
fortunately is definitely not the case. There is clear evidence that instrumental 
biases from the satellite observations and the redshift measurements result in 
a marked distortion of the 'true' GRB redshift distribution. 

Only about 50% of Swift GRBs have been observed with an optical afterglow, 
despite the rapid localization and follow up by telescopes capable of spec- 
troscopy. This leaves a significant fraction of the observed population with 
no possibility of observing a redshift directly from the afterglow. One possi- 
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ble reason for this deficiency is dust extinction of the afterglow i n the GRB 



host g alaxy, but this may not account for all of the deficiencies. 



Fiore et al 



( 120061 ) showed that the different redshift distributions obtained using Bep- 
poSax/HETE2 and Swift were at least partially due to the sensitivity of the 
different instruments, localisation efficiency and efficiency of obtaining a spec- 
troscopic redshift of the afterglow or of the host galaxy. The afterglow bright- 
ness is an important factor in obtaining a spectroscopic redshift. The redshifts 
that have been measured come from bright GRBs, implying the measured 
GRB redshift distribution is a sub-population of the bright part of the GRB 
LF. 



An important question is how selection biases affect the fitting procedures 
employed to extract the SFR and/or GRB rate evolution using redshift data? 
If this bias is not accounted for in the GRB scaling f unction , the m odel for 



the observed redshift distribution will be biased. The [Band 



( 120061 ) analysis 



showed that Swift has increased sensitivity to high redshift bursts because 
of their longer duration. This selection effect may manifest in the redshift 
distribution as an evolution of the GRB rate as a function of redshift. The 
inability to account for these biases in the model fits to the redshift distribution 
could easily result in misinterpreting fitting results. In addition, the high-z 
regime where this may be a problem is also where there is the most uncertainty 
in the GRB rate density and SFR. 



To account for the most significant biases requires a 'complete GRB selection 
function'. An expression for such a function could consist of a product of 
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normalized scaling functions for the main classes of biase 

^complete (^) = ^GRB (^satellite (^redshift (^) , (13) 

where ^satellite (z) would account for the bias discussed by Band (2006) and 
V'redshift^) would incorporate a model for the 'redshift dropout' from the dif- 
ficulty of obtaining resolved lines in different redshift regimes (Firore et al. 
2006). Various attempts have been made to at least minimize these biases, 
including selecting only those redshifts that pass certain criteria, such as the 
source position in relation to Galactic dust obscuration. 

The use of Monte Carlo simulations to study the propagation of uncertainty 
from the GRB selection function to the fitting parameters for GRB evolution 
will provide a better understanding of the uncertainty in the best-fit models. 
This has already been performed using different LFs, but has not been fully 
extended to instrumental biases, partly because of the complexity of the biases. 
Nonetheless, unless the complete GRB selection function can be constrained to 
a better precision, it is unlikely that it will be possible to discriminate between 
different GRB evolutionary scenarios based on the redshift distribution alone. 

5 Concluding remarks 

Presently there is no 'complete' GRB selection function to account for the bi- 
ases emerging in the GRB redshift data. Nonetheless there is an urgent need to 
quantify selection effects that plague the observed GRB redshift distribution. 

2 We note the complete scaling function, Complete (z), or its components, iPgrb(z) 
and ^satellite (z), are sometimes termed 'detection efficiency' curves. 
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The result of this uncertainty in the redshift distribution is that it is difficult 
to differentiate between an evolving GRB luminosity function, enhanced rate 
density and the shape of the SFR in z > 2 — 3. It is very difficult to place 
quantitative confidence limits on the fitted GRB evolution parameters using 
the redshift distribution because the uncertainties inherent in the models are 
not fully quantified, especially VWshift^)- Previous studies exploit both the 
brightness and redshift distributions to explore different GRB evolutionary 
scenarios and estimate parameter limits based on a range of plausible mod- 
els. If all the instrumental biases from the redshift distribution are accounted 
for, we can explore the physical GRB rate evolution parameter space with 
confidence and our main limitation will be small number statistics, a problem 
that will inevitably be solved by acquiring more data from Swift and other 
satellites. 
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3 4 5 6 7 




0.4 
log(l+z) 



Fig. 1. Evolution of the SFR density with redshift reproduced from 



Hopkins fc Beacoml (|2006l ). They have compiled data from various sources (see ref- 
erences therein) including FIR (24 jiim) — hatched region & triangles, UV — squares, 
radio(1.4 GHz) — open red star and UDF estimates — crosses. The solid lines are 



their fits to these data using parametrized fitting functions. 
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Fig. 2. Plots of the cumulative redshift distributions using 42 and 48 pre-Swift and 
Swift bursts respectively. The median redshift of the two distributions are 1.0 and 
2.4 respectively, indicating the presence of different biasses. 
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redshift z 

Fig. 3. A plot of the GRB scaling function, iPgrb( z ), using the LF expressed in 
equation ([2]) and the peak observable flux from equation (J3]) assuming a limiting 
detector sensitivity similar to Swift of 10 -8 cm~ 2 s _1 . The Band function with 
spectral slopes a = — 1 and j3 = —2.25, integrated over the 50-300 keV band, is 
used to define the source luminosity. 
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8 



Fig. 4. Plot of the probability distribution of GRBs as a function of redshift, assum- 



ing the scaling function plotted in figure El GRB rate evo 
the SFR and for definiteness we use a SFR model from 



ution is assumed to 



:bllow 



Porciani fc Madaul ((200l|) 



that peaks in z = 1 — 2, and remains approximately constant at z > 2. The model 
curves do not account for selection biases associated with the satellite detectors or 
the redshift determinations. 



26 



